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guide change  markedly  over  lateral  distances  comparable  to  a wavelength  or  Fresnel 
zone.  Full-wave  solutions  are  used  to  describe  both  the  vertical  and  horizontal 
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describing  the  vertical  dependence, l this  report  focuses  on  solving  the 
equation  for  the  lateral  dependence. 'A To  facilitate  numerical  solution,  the 
lateral  equation  is  transformed  into  an  integral  equation  that  accounts  for 
most  full-wave  properties,  including  diffraction  around  a localized  disturbance 
and  reflection  from  lateral  gradients.  Numerical  solutions  based  on  model 
disturbances  having  lateral  gradients  in  the  direction  of  propagation  reveal 
a standing  wave  pattern  in  front  of  the  disturbance.  The  pattern  is  pronounced 
if  the  waveguide  properties  change  substantially  ove&a  distance  equal  to 
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The  often  used  WKB  approximation  omits  the  standing1 wave  pattern  and 
thus  gives  poor  accuracy  for  regions  in  front  of  a disturbance  having  a 
relatively  abrupt  boundary.  Because  it  ignores  gradient  reflections,  the 
WKB  approximation  also  overestimates  the  signal  transmitted  beyond  a non- 
uniformity in  the  waveguide.  However,  for  all  models  considered — abrupt 
or  diffuse — the  overestimate  is  minor. 


HWriASSIFIED 

L ASSIFlCATlON  OF 


••CUeiTY  CLASSIFICATION  of  this  »AOf'n»n  Of  *•>«•«**) 


11 


PREFACE 

This  report  develops  a method  for  analyzing  the  propagation  of  ELF 
(extremely  low-frequency)  waves  in  the  presence  of  laterally  nonuniform 
ionospheric  disturbances  that  violate  the  validity  criteria  of  the 
widely  used  WKB  approximation.  The  method  is  particularly  useful  for 
calculating  ELF  propagation  anomalies  caused  by  nuclear  detonations  at 
altitudes  of  30  to  150  km. 


| The  present  report  continues  the  Pacific-Sierra  Research  Corporation's 

analysis  of  long-wave  propagation  in  nuclear  and  naturally  disturbed 
environments.  In  focusing  on  the  effects  of  lateral  ionospheric  gradients 
^ in  the  direction  of  propagation,  it  complements  an  earlier  investigation 

of  the  effects  of  gradients  transverse  to  the  propagation  path  (E.  C.  Field, 


* ELF  Propagation  in  a Non-Strati fied  Earth-Ionosphere  Waveguide , PSR  806, 

| April  1978). 


iii 


! 


I 


l 


» 


I 3 


C 


ft 


3 


SUMMARY 

This  report  develops  and  applies  a practical  method  for  analyzing  long- 
wave propagation  under  conditions  where  the  properties  of  the  earth- ionosphere 
waveguide  change  markedly  over  lateral  distances  comparable  to  a wavelength 
or  Fresnel  zone.  Full-wave  solutions  are  used  to  describe  both  the  vertical 
and  horizontal  dependences  of  the  fields,  but  certain  compromises  are  made 
to  achieve  tractability.  The  method  is  thus  characterized  as  "quasl-full- 
wave."  Its  main  limitation  is  that  the  equation  describing  the  vertical 
dependence  is  assumed  to  nearly  decouple  from  the  equation  governing  the 
lateral  dependence.  The  method  is  valid  at  any  frequency  for  which  waveguide 
modes  describe  terrestrial  propagation.  Nonetheless,  its  practical  utility 
is  probably  limited  to  ELF  because  the  lateral  properties  of  the  earth- 
ionosphere  waveguide  are  usually  gradual  enough  to  permit  use  of  the  WKB 
approximation  at  higher  frequencies. 

Since  the  literature  provides  abundant  solutions  for  the  equation 
describing  the  vertical  dependence,  this  report  focuses  on  solving  the 
equation  for  the  lateral  dependence.  To  facilitate  numerical  solution, 
the  lateral  equation  is  transformed  into  an  integral  equation  that  accounts 
for  most  full-wave  properties,  including  diffraction  around  a localized 
disturbance  and  reflection  from  lateral  gradients.  Numerical  solutions 
based  on  model  disturbances  having  lateral  gradients  in  the  direction  of 
propagation  reveal  a standing  wave  pattern  in  front  of  the  disturbance. 

The  pattern  is  pronounced  if  the  waveguide  properties  change  substantially 
over  a distance  equal  to  about  one-sixth  of  a wavelength,  and  is  minor  if 
the  disturbance  is  more  diffuse. 

The  often  used  WKB  approximation  omits  the  standing  wave  pattern  and 


thus  gives  poor  accuracy  for  regions  in  front  of  a disturbance  having  a 
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relatively  abrupt  boundary.  Because  it  ignores  gradient  reflections,  the 
WKB  approximation  also  overestimates  the  signal  transmitted  beyond  a 
nonuniformity  in  the  waveguide.  However,  for  all  models  considered — 
abrupt  or  diffuse — the  overestimate  is  minor. 
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I.  INTRODUCTION 


The  wavelength  and  first  Fresnel  zone  for  ELF  waves  propagating  in  the 
earth- ionosphere  waveguide  can  extend  several  megameters.  Many  types  of 
ionospheric  irregularities  exhibit  significant  lateral  variations  over  such 
a distance.  Ideally,  the  effects  of  such  irregularities  on  ELF  propagation 
should  be  analyzed  by  full-wave  methods  rather  than  the  often  used  WKB  approx- 
imations, which  work  well  for  higher  frequency  signals  having  shorter  wave- 
lengths. To  our  knowledge,  however,  complete  full-wave  solutions  that 
simultaneously  account  for  both  vertical  and  lateral  ionospheric  inhomo- 
geneitles  have  to  date  proven  intractable. 

This  report  develops  and  applies  a practical  method  for  analyzing  long- 
wave propagation  under  conditions  where  the  properties  of  the  earth-ionosphere 
waveguide  change  markedly  over  lateral  distances  comparable  to  a wavelength 
or  Fresnel  zone.  Full-wave  solutions  are  used  to  describe  both  the  vertical 
and  horizontal  dependences  of  the  fields,  but  certain  compromises  are  made 
to  achieve  tractabllity.  The  method  is  thus  characterized  as  "quasi-full- 
wave."  Its  main  limitation  is  that  the  equation  describing  the  vertical 
dependence  is  assumed  to  nearly  decouple  from  the  equation  governing  the 
lateral  dependence.  That  assumption  appears  reasonable  for  ionospheric 
irregularities  of  large  enough  extent  to  significantly  affect  long-wave 
propagation.  The  vertical  equation  can  be  solved  using  the  well-known 
method  developed  by  Budden  (1961a)  and  often  applied  in  the  literature. 

This  report  focuses  on  solving  the  equation  for  the  lateral  dependence  of 
fields.  To  facilitate  numerical  solution,  the  lateral  equation  is  trans- 
formed into  an  Integral  equation  analogous  to  that  derived  by  Hufford  (1962) 
to  describe  ground-wave  propagation  over  irregular  terrain. 
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A previous  report  (Field,  1978)  developed  approximate  solutions  for 
a weak,  localized  ionospheric  disturbance  re- ote  from  the  transmitter  and 
receiver.  It  pointed  up  the  importance  of  gradients  transverse  to  the 
propagation  path.  The  present  report  complements  the  earlier  results  by 
giving  solutions  for  ionospheric  disturbances  of  arbitrary  strength  and 
extent  that  are  azimuthally  symmetric  about  the  transmitter. 

Section  II  derives  the  Integral  equation  that  describes  the  lateral 
dependence  of  the  fields.  Section  III  gives  numerical  results  for  propagation 
in  the  presence  of  ionospheric  disturbances  centered  over  the  transmitter  and 
over  the  midpath  of  a long  ELF  link.  Several  of  the  model  disturbances 
represent  environments  produced  by  single,  high-altitude  nuclear  bursts. 

All  results  are  compared  with  WKB  solutions,  and  the  effects  of  varying 
the  lateral  scale  size  and  boundary  diffuseness  of  the  assumed  disturbances 
are  examined.  The  conclusions  are  presented  in  Sec.  IV. 
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II.  INTEGRAL  EQUATION  FOR  LATERAL  DEPENDENCE  OF  FIELDS 
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First  we  review  the  well-known  equations  for  ELF  propagation  under 
stratified  conditions,  then  derive  the  integral  equation  that  describes 
the  lateral  dependence  of  the  fields  for  nonstratif led  conditions.  Two 
simpler  versions  of  the  integral  equation  are  derived  for  the  special 
cases  of  a localized  disturbance  remote  from  the  terminals  and  an  azimu- 
thally  symmetric  disturbance  centered  over  the  transmitter.  To  avoid 
mathematical  complexities  unrelated  to  lateral  ionospheric  gradients,  we 
ignore  earth  curvature  and  the  geomagnetic  field.  The  first  approxi- 
mation is  well  justified  at  ELF;  the  second  yields  results  that  are 
reasonably  accurate  for  ambient  daytime  conditions  and  very  accurate  for 
disturbances  where  ionospheric  reflection  heights  are  depressed  below 
ambient  levels.  The  equations  derived  below  are  inappropriate  for  esti- 
mating propagation  under  normal  nighttime  conditions.  A time  dependence 
e*wt  is  assumed. 

STRATIFIED  CONDITIONS 

We  begin  by  defining  a function,  <|i,  that  describes  the  lateral 
dependence  of  the  fields  in  the  earth-ionosphere  waveguide.  For  ELF 
TEM-mode  propagation  at  ranges  exceeding  the  effective  ionospheric 
reflection  height,  <| > could  denote  either  the  z component  of  the  electric 
field  or  vector  potential,  or  the  horizontal  component  of  the  magnetic 
intensity.  Here  we  associate  <(i  with  the  vertical  electric-field 
component,  E. 

ELF  signals  are  typically  radiated  from  a horizontal  electric 
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dipole  antenna.  Hence  the  fields  have  both  radial  and  azimuthal 
dependences,  even  if  the  earth-ionosphere  waveguide  is  laterally 


» 

» 


homogeneous.  Under  typical  conditions,  then,  the  wave  equation  can  be  solved 
by  separation  of  variables,  which  gives  the  following  well-known  result  (e.g.. 


Galejs,  1972): 


Eq  - AAoFQ(z)i(;o(r,0)  volts/m 


In  Eq.  (1)  A is  a constant  involving  dipole  moment,  wave  frequency,  and 
ground  conductivity;  Aq  is  the  excitation  factor  describing  the  efficiency 
with  which  the  TEM  mode  is  launched;  and  Fq(z)  is  the  height-gain  function 
describing  the  vertical  dependence  of  the  field,  normalized  to  unity  at 
z ■ 0.  Throughout  this  report  the  subscript  0 is  used  to  denote  quantities 
associated  with  undisturbed,  laterally  homogeneous  conditions. 

The  function  <||q  in  Eq.  (1)  satisfies  the  two-dimensional  wave  equation 

(V*  + k2SoH0  * o I 


where  VT  is  the  transverse  Laplacian,  k is  the  free-space  wave  number,  and 
Sq  is  a propagation  constant  determined  by  imposing  boundary  conditions  on 
Fq  at  the  ground  and  in  the  ionosphere.  For  stratified  conditions,  Eq.  (2) 
is  easily  solved  to  give 


4,q  - cos0  SqHj  (kSQr) 


where  r and  9 are  the  usual  circular  cylindrical  coordinates;  the  horizontal 
dipole  transmitter  is  located  at  r * 0 and  oriented  at  0 * 0;  and  H is  the 
Hankel  function. 


NONSTRATIFIED  CONDITIONS 

In  the  presence  of  lateral  ionospheric  gradients  a rigorous  separation 
of  variables  is  impossible.  Hence,  the  fields  cannot  be  expressed  as  a 
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product  of  vertical,  azimuthal,  and  radial  functions  as  in  Eqs.  (1)  and  (3). 
However,  for  ionospheric  inhomogeneltles  large  enough  to  significantly 
affect  ELF  propagation,  scale  lengths  for  lateral  variations  of  the 
ionospheric  refractive  index  tend  to  be  at  least  an  order  of  magnitude 
larger  than  those  for  vertical  variations.  It  can  therefore  be  argued 
on  semiquantitatlve  grounds  that  the  vertical  dependence,  F,  of  the  fields 
and  the  associated  eigenvalue,  S,  are  governed  primarily  by  the  local 
ionosphere.  In  this  approximation  it  is  assumed  by  analogy  with  Eq.  (1) 
that  the  ground-level  field  is  given  by 

E ss  AA(x,y)«|»(x,y)  . (4) 

To  account  for  lateral  variations  that  are  nonsymmetric  about  the  source, 
we  have  switched  from  cylindrical  to  Cartesian  coordinates.  Equation  (4) 
also  reflects  the  fact  that  F * 1 at  z * 0. 

To  find  the  field  from  Eq.  (4),  we  use  a mixture  of  eikonal  and  full- 
wave  techniques  perhaps  best  described  as  a "quasi-full-wave"  method.  By  the 
arguments  above,  the  ionosphere  is  assumed  locally  stratified  for  finding  the 
vertical  dependence  of  the  fields,  the  eigenvalue  S(x,y),  and  the  excitation 
factor,  A.  Thus  the  equation  for  these  quantities,  and  the  method  of  solution, 
are  formally  Identical  to  those  widely  applied  in  analyzing  ELF  propagation  in 
a laterally  uniform  earth- ionosphere  waveguide.  They  differ  in  practice, 
however,  because  the  equation  for  the  vertical  dependence  must  be  solved  at 
a large  number  of  locations,  (x,y),  each  characterized  by  a local  ionospheric 
height-profile.  On  the  other  hand,  a single  solution  suffices  for  all 
locations  under  laterally  uniform  conditions. 

The  lateral  dependence  of  the  field  is  expressed  mainly  in  the  function 
<»,  which  by  analogy  with  Eq.  (2)  satisfies  the  equation 


(V2  + k2S2  (x,y)  )4>  - o 


The  local  eigenvalue  S(x,y),  unlike  the  propagation  constant  S^,  exhibits 

spatial  dependence,  which  precludes  general  analytic  solutions  of  Eq.  (5). 

S(x,y)  is  found  by  imposing  boundary  conditions  on  the  fields  in  the 

ionosphere  and  at  the  ground  for  a large  number  of  geographic  locations, 

thus  obtaining  a matrix  of  values  to  be  entered  in  Eq.  (5).  The  literature 

supplies  full-wave  methods  for  obtaining  S and  A for  virtually  any  ionospheric 

re tractive- index  height -profile,  as  well  as  numerical  results  for  many 

ambient  and  disturbed  models  of  the  ionosphere  (e.g. , Budden , 1961a ; Field, 

1970;  Wait , 1970;  Galejs,  1972;  Pappert  and  Moler,  1974;  Grei'finger  and 

Grei finger,  1978).  We  can  therefore  assume  that  all  quantities  except  i|> 

are  either  known  or  readily  obtainable,  which  allows  us  to  concentrate  on 

obtaining  full-wave  solutions  of  Eq.  (5)  for  the  lateral  wave  function,  if/. 

To  facilitate  numerical  solution,  we  first  recast  the  problem  into 

an  integral  equation.  Our  main  interest  is  in  calculating  the  effects  on 

ELF  propagation  of  a laterally  nonuniform  ionospheric  disturbance,  which 

2 2 

can  be  characterized  by  the  difference  S (x,y)  - S^.  The  undisturbed  wave 
function,  <J»q,  is  governed  by  Sq  and  can  be  assumed  known  through  Eq.  (3). 
Following  Wait  (1964),  the  subtraction  of  Eq.  (2)  from  Eq.  (5)  leads  to 


(v2  + k2sj)  (♦-*(,)  - -k2[s2-sj]* 


We  convert  Eq.  (6)  to  the  desired  integral  equation  by  using  the  Green's 
function. 


0 - -l.«<!)(kS0r2> 


(7) 
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which  satisfies 


(V2  + k2Sp)G  - -4n6(r2) 


(8) 


where  r2  ■ | r - r j | (r  and  are  vectors  from  the  origin  to  the  observa- 
tion point  and  to  an  integration  point,  respectively).  By  following  the  usual 
Green '8  function  procedure  and  applying  the  two-dimensional  Green's  theorem, 
the  Integral  equation  for  is  obtained: 

«|»(x,y)  - i|»0 (x,y)  ~ f JA*' dy' fS2(x’ ,y')  - S2  ]Hq2)  (IcSq^iKx’  ,y' ) . (9) 

Equation  (9)  is  formally  identical  to  an  integral  equation  given  by  Uait 
(.1964),  who  used  first-order  perturbation  theory  to  obtain  approximate 
solutions  valid  at  VLF  (very  low  frequencies). 

It  is  also  convenient  to  define  a relative  propagation  function,  U, 
which  denotes  the  fractional  amount  by  which  the  disturbed  lateral  wave 
function,  differs  from  the  undisturbed  function,  t|i^.  Specifically,  we 
define  the  relative  propagation  function  by 


♦ (x,y)  = W(x,y)i|»0(x,y) 


(10) 


The  insertion  of  Eqs.  (10)  and  (3)  into  Eq.  (9)  gives  the  following  integral 
equation  for  W: 


2 * 

W(x,y)  - 1 - yy,dx'dy'[s2(x',y')  - sj] 


r«^l  H02)<kV2>  ■l2)<“0r1)  w(x,  .. 

W H<2>(kS0r)  ,y  ’ 


(ID 
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where  r2  - x2  + y2,  r2  - (x')2  + (y')2,  and  r2  ■ (x  - x')2  + (y  - y')2. 

Equation  (11)  Is  the  most  general  form  of  the  Integral  equation  for  the 

relative  propagation  function.  Note  that  in  the  absence  of  an  ionospheric 
2 2 

disturbance,  S * Sq  and  Eq.  (11)  has  the  trivial  solution  W « 1. 

LOCALIZED  DISTURBANCES 

Despite  the  great  wavelengths  at  ELF,  the  distance,  r,  between  trans- 
mitter and  receiver  is  usually  large  enough  to  permit  use  of  the  asymptotic 
approximation  for  the  Hankel  function  in  the  denominator  of  Eq.  (11).  In 
the  numerator,  however,  the  validity  condition  on  the  asymptotic  approxima- 
tion for  the  Hankel  functions  is  much  more  restrictive.  That  is 

because  r^  or  r2  can  be  very  small  over  part  of  the  integration  interval 
2 2 

if  S - Sq  is  nonzero  near  the  transmitter  or  receiver.  Thus,  we  may  use 

the  asymptotic  formula  only  if  the  disturbance  occurs  at  least  a megameter 

2 2 

away  from  either  terminal;  S - Sq  cannot  be  nonzero  within  a megameter  of 
the  transmitter  or  receiver.  If  that  condition  is  satisfied,  the  following 
asymptotic  approximation  of  Eq.  (11)  may  be  used: 


W(x,y)  - 1 + 


fc3'2.-1'* 

2/2ls: 


00 

ff  d*'dy'[s2  - s2] 


[^[37]  “!>|-lkso(v  r2  - ']!  • 


Not  surprisingly,  Eq.  (12),  which  describes  propagation  in  a laterally 
Irregular  earth-ionosphere  waveguide,  is  formally  similar  to  the  classic 
integral  equation  representing  ground-wave  propagation  over  Irregular 
terrain  (e.g. , Huffotvl,  1952). 
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DISTURBANCE  SYMMETRIC  ABOUT  THE  TRANSMITTER 

If  an  azlmuthally  symmetric  disturbance  is  centered  over  the  trans- 
2 2 

mitter,  the  quantity  S - SQ  in  Eq.  (9)  depends  only  on  the  distance. 


2 2 

(x')  + (y')  , from  the  origin.  For  that  situation  Eq.  (9)  can  be  rewritten 


in  cylindrical  coordinates,  and  the  azimuthal  dependence  can  be  Integrated 
out  to  give  a one-dimensional  integral  equation  for  the  relative  propagation 
function.  (The  procedure  is  outlined  in  Appendix  A.)  The  resulting  integral 
equation  is 


I C 


r. 


W(r)  - 1 - 


ik2x 


OD 

/ 


dr'K(r,r')W(r’) 


r'  < r 


K(r,r*)  - 


r’[S‘(r’)  - S2]  J j (kSQr)Hj2) (kSQr ' )H^2) (kSQr ' ) 


H^kS^) 


r'  > r 


J denotes  the  Bessel  function  in  Eq.  (13b). 

In  the  next  section,  the  full-wave  results  computed  from  Eq.  (13) 
are  compared  with  the  often  used  WKB  approximation,  which  is  valid  at 
large  distances  from  the  source  for  lateral  inhomogeneities  with  scales 
greater  than  an  Inverse  wave  number.  As  shown  in  Appendix  B,  the  WKB 
approximation  of  the  relative  propagation  function  for  a disturbance 
synsetric  about  the  transmitter  is  given  by 


W(r)  ^ 
WKB 


exp  ~lk 

/ (S(r')  - S_ J dr ' 

y(r>J  l o 

f ^ 

(13a) 


where  r is  the  distance  from  the  transmitter  to  the  observation  point  and 
r'[S2(r’)  - S2  ] J j (kSQr ' )H^2> (kSQr ' ) 


(13b) 


(14) 
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III.  NUMERICAL  RESULTS 

This  section  presents  full-wave  numerical  solutions  of  the  one- 
dimenslonal  Eq.  (13),  which  applies  to  ionospheric  disturbances  that 
depend  only  on  distance  from  the  transmitter.  The  solutions  illustrate 
the  effects  of  gradients  in  the  direction  of  propagation.  We  have  not 
yet  developed  numerical  procedures  for  solving  the  two-dimensional  Eq.  (11), 
which  applies  to  disturbances  that  depend  on  both  lateral  coordinates. 
However,  as  mentioned  in  Sec.  I,  Field  (1978)  gives  approximate  solutions 
of  Eq.  (12)  for  weak,  localized  disturbances  that  depend  on  two  lateral 
coordinates.  Those  solutions  pertain  to  gradients  transverse  to  the 
direction  of  propagation  and  thus  complement  the  solutions  given 
below. 

To  check  on  the  accuracy  of  the  algorithm  used  to  solve  Eq.  (13),  we 
performed  numerical  calculations  for  two  idealized  types  of  disturbance 
for  which  (and  hence  W)  could  be  analytically  found  from  Eq.  (5).  The 
idealized  types  are  l)  a laterally  uniform  disturbance  where  S is  indepen- 
dent of  r but  differs  from  S^,  and  2)  a nonuniform  disturbance  where 
2 2 2 

S - Sq  varies  as  1/r  . The  solution  to  Eq.  (5)  is  a simple  Hankel  func- 
tion for  the  first  test  case  and  a complicated  combination  of  Hankel 
functions  of  complex  order  for  the  second.  Both  solutions  satisfy  Eq.  (13) 
and  agree  almost  exactly  with  the  numerical  solutions  of  Eq.  (13). 

MODELS  USED  FOR  NUMERICAL  CALCULATIONS 

To  illustrate  the  dependence  of  the  relative  propagation  function  on 
the  severity  of  gradients  in  the  direction  of  propagation,  we  use  three 
simple  but  realistic  models  for  S.  Figure  1 depicts  the  models  schematically 
and  gives  the  corresponding  analytic  expressions.  In  each  model  the  distur- 
bance is  characterized  by  an  exponential  function  with  a scale  length,  Ar, 


! 


11 


r(Mm) 

a:  Disturbance  over  transmitter  ((S2-Sq)  * (S2-Sq)  exp  (-r2/(Ar)2)) 


r|Mm) 


b:  Disturbance  over  midpath  ((S2-Sq)  = (S^-Sq)  exp  (-(r-5)2/(Ar)2>) 


c:  Uniform  regions  separated  by  diffuse  boundary  at  midpath 
((S-Sq)  - (Sj-Sgl/d+expl-fr-SJ/Ar])) 

Fig.  1— Models  of  azlmuthally  symmetric  disturbances  used  In 
calculations  (not  to  scale). 


that  can  be  varied  to  indicate  either  abrupt  or  diffuse  boundaries.  The 
first  model  (Fig.  la)  represents  an  ionospheric  depression  centered  directly 
over  the  transmitter.  The  second  model  (Fig.  lb)  represents  a depression 
centered  over  the  midpath  of  a 10-Mm  link.  The  third  model  (Fig.  lc) 
represents  propagation  from  a uniformly  undisturbed  region  into  a region 
of  uniformly  depressed  ionosphere,  the  two  regions  being  separated  by  a 
diffuse  boundary  centered  at  midpath.  Because  of  the  large  radius  of 
curvature,  the  models  in  Figs,  lb  and  lc  should  reasonably  represent  a 
plane  waveguide-mode  incident  on  a one-dimensional  disturbance.  Note  that 
a lowering  of  the  "effective"  ionospheric  height  usually  increases  both 
the  real  and  imaginary  parts  of  S. 

The  symbols  in  Fig.  1 are  defined  as  follows:  SQ  is  the  eigenvalue 
governing  propagation  in  the  undisturbed  region;  Sj  is  the  eigenvalue 
governing  propagation  at  the  most  disturbed  point;  and  Ar  (Mm)  is  the  charac- 
teristic length  over  which  the  conditions  change  from  undisturbed  to  disturbed, 
i.e.,  S(r)  effectively  changes  from  Sq  to  over  a distance  two  or  three 

times  as  large  as  Ar.  The  table  below  shows  the  numerical  values  and  corres- 
ponding attenuation  rates  (in  dB/Mm)  used  for  the  constants  and  Sj.  Sq 
values  represent  nominal  ambient  daytime  conditions;  those  for  Sj  represent 
a severe  ionospheric  disturbance,  such  as  a major  solar  proton  event. 


Table 

PARAMETERS  USED  IN  CALCULATIONS 


f(Ha) 

Sq  (undisturbed) 

Sj  (disturbed) 

45 

1.2  - 0.08  i 

1.75  - 0.3  i 

(0.65  dB/Mm) 

(2.4  dB/>ta) 

75 

1.15  - 0.085  i 

1.5  - 0.25  i 

(1.15  dB/no 

(3.4  dB/Mm) 
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The  simple  analytic  forms  In  Fig.  1 are,  of  course,  used  for  Interpre- 
tive convenience.  The  numerical  methods  employed  to  obtain  the  corresponding 
solutions  to  Eq.  (13)  work  equally  well  for  any  disturbance  where  S depends 
only  on  r.  Then,  however,  S(r)  must  be  calculated  numerically  and  entered 
as  a table  to  the  algorithm  for  solving  Eq.  (13).  To  illustrate  this 
capability,  results  are  also  given  for  ELF  propagation  in  the  presence  of 
a disturbance  created  by  detonation  of  a high-altitude  nuclear  weapon 
directly  over  the  transmitter. 

CALCULATED  RESULTS 

All  results  given  below  pertain  to  the  relative  propagation  function, 

W,  the  ratio  of  the  disturbed  to  undisturbed  radial  wave  function,  <|t.  Recall 
that  the  ground-level  electric  field  is  proportional  to  the  product  of  and 
the  excitation  factor,  A (see  Eq.  (A)).  W does  not  contain  the  excitation 
factor  and  therefore  does  not  always  represent  the  ratio  of  the  disturbed  to 
undisturbed  fields.  The  excitation  factor  is  roughly  proportional  to  the  in- 
verse  of  the  thickness  (i.e.,  effective  height)  of  the  earth-ionosphere  wave- 
guide. For  laterally  nonuniform  waveguides,  the  geometric  mean  of  the  Inverse 
thickness  at  the  transmitter  and  receiver  is  often  used  as  an  approximation. 
Ionospheric  disturbances  usually  depress  the  ionosphere,  thereby  Increasing  the 
local  excitation  factor.  Thus,  if  either  the  transmitter  or  receiver  is  in  a 
disturbed  region,  as  in  Figs,  la  and  lc,  the  excitation  factor  is  larger 
than  the  ambient  value  and  the  ratio  of  disturbed  to  undisturbed  fields  is 
larger  than  the  ratio  W.  However,  W does  represent  the  ratio  of  fields  if 
both  terminals  are  In  undisturbed  regions,  as  depicted  in  Fig.  lb.  Regarded 
heurlstically,  W accounts  for  changes  in  attenuation  and  phase  velocity, 
whereas  A accounts  for  the  fact  that  the  power  density  in  the  wavefront  is 
Inversely  proportional  to  the  cross-sectional  area  of  the  waveguide. 
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The  calculated  results  for  the  ionospheric  depression  centered  over 

the  transmitter,  illustrated  in  Fig.  la,  are  shown  in  Figs.  2 through  4. 

Figure  2 shows  the  magnitude  of  W versus  the  lateral  scale,  Ar,  of  the 

depression  for  a pathlength  of  10  Mm  and  frequencies  of  45  and  75  Hz. 

Curves  depict  the  full-wave  solutions  to  the  integral  equation  (Eq.  13) 

and  the  approximate  WKB  solution  given  by  Eq.  (14).  As  expected,  the 

full-wave  and  WKB  solutions  agree  closely  if  Ar  > X/2it,  where  X is  the 

* 

free-space  wavelength.  This  condition  satisfies  the  familiar  WKB  validity 
criterion  that  the  waveguide  properties  change  only  slightly  over  a hori- 
zontal distance  equal  to  an  inverse  wave  number.  For  smaller  values  of 

1 

Ar  that  correspond  to  highly  localized  disturbances  having  relatively  large 
horizontal  gradients,  the  WKB  solutions  significantly  overstate  the  magnitude 
of  W. 

Figure  3,  analogous  to  Fig.  2,  shows  the  phase  of  W as  a function  of 
Ar  for  a disturbance  over  the  transmitter.  Again,  agreement  between  the 
full-wave  and  WKB  solutions  is  8°od  provided  the  horizontal  scale  of  the 
disturbance  exceeds  a megameter  or  so. 

Figure  4 shows  the  magnitude  of  W as  a function  of  distance  rather 
than  scale  size,  Ar,  as  in  Fig.  2.  Recall  that  W is  the  ratio  of  the 
disturbed  to  undisturbed  radial  wave  functions,  and  increases  or  decreases 
in  electric-field  strength  are  characterized  by  W's  greater  or  less  than 
unity,  respectively.  Figure  4 shows  that  widespread  disturbances  (large  Ar) 
cause  a net  reduction  in  the  radial  wave  function,  whereas  more  confined 
disturbances  actually  increase  it  slightly.  (Of  course,  as  Ar  approaches 
zero  the  disturbance  essentially  ceases  and  W approaches  unity.) 


© 


The  quantity  X/2x  is  about  1 Ma  at  45  Hz  and  about  0.6  Mm  at  75  Hz. 


f W versus  distance  for  disturbance  over  transmitter 
a and  table). 
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This  phenomenon  is  a result  of  two  competing  factors:  1)  increased 
attenuation,  which  causes  U to  decrease  with  distance,  and  2)  propagation 
from  a disturbed  region  of  low-phase  velocity  (large  S)  into  an  undisturbed 
region  of  higher  phase  velocity  (small  S) , which  causes  the  electric  field, 
and  hence  W,  to  increase  with  distance.  The  latter  is  analogous  to  the 
well-known  increase  of  the  electric  field  that  occurs  with  a reduction  in 
the  refractive  index  of  a propagation  medium.  Clearly,  the  importance  of 
anomalous  attenuation  is  proportional  to  the  pathlength  exposed  to  the  dis- 
turbance, whereas  any  increase  in  the  electric  field  depends  mainly  on 
the  contrast  between  phase  velocities  at  the  transmitter  and  receiver.  Thus, 

W is  slightly  increased  by  localized  disturbances  with  small  exposed  path- 
lengths,  whereas  the  degrading  effects  of  increased  attenuation  dominate  in 
widespread  disturbances.  The  curve  labeled  Ar  = 00  in  Fig.  4 shows  how  W is 

. j 

affected  by  a horizontally  uniform  disturbance  (the  value  Sj  in  the  table) 
at  all  ranges.  For  this  limiting  case,  W decreases  essentially  monotonically 
because  the  phase  velocity  is  independent  of  distance,  whereas  the  attenuation 
rate  everywhere  increases  over  the  ambient  value. 

Figures  5 through  7 show  the  calculated  values  of  W versus  distance  for 
the  disturbance  at  r **  5 Mm  illustrated  in  Fig.  lb.  These  figures  portray  the 
effects  of  making  the  disturbance  progressively  broader  and  hence  less  abrupt. 

A distinct  standing  wave  pattern  due  to  reflections  from  the  disturbance  is 
evident  in  Figs.  5 and  6,  which  show  results  for  the  two  cases  (Ar  ■ 0.2  and 
0.5  Mm)  where  the  waveguide  properties  change  substantially  over  a distance 
of  X/2x.  Since  the  WKB  solution  ignores  reflections,  it  incorrectly  omits 
the  standing  wave  pattern  in  front  of  the  disturbance  and  overestimates  the 
signal  transmitted  through  the  disturbance.  Nevertheless,  even  for  distur- 
bances as  abrupt  as  those  in  Figs.  5 and  6,  the  WKB  solution  gives  a remark- 
ably good  approximation  of  the  signal  behind  the  disturbance.  In  Fig.  7 


--Amplitude  of  W versus  distance  for  disturbance  at 


-Amplitude  of  W versus  distance  for  disturbance  at  r s 5 Wn: 
Ar  = 0.5  Hn,  f = 75  Hz  (see  Fig.  lb  and  table). 


Amplitude  of  W versus  distance  for  disturbance  at 
Ar  = 1 Hn,  f = 75  Hz  (see  Fig.  lb  and  table). 
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the  disturbance  (Ar  ■ 1 Mm)  changes  fairly  slowly  over  the  distance  X/2n. 

As  expected,  the  standing  wave  pattern  is  greatly  diminished,  and  the 
WKB  solution  affords  an  excellent  approximation  at  all  ranges. 

Figure  8 illustrates  propagation  from  a uniformly  undisturbed  region 
into  a uniformly  disturbed  one  (see  Fig.  lc)  where  the  boundary  between  the 
regions  is  quite  abrupt  (Ar  - 0.2  Mm)  relative  to  the  distance  X/2n.  Again,  a 
standing  wave  pattern — absent  from  the  WKB  solution — is  evident  in  front 
of  the  boundary,  and  W falls  off  monotonically  behind  the  boundary  because 
of  the  increased  attenuation  associated  with  the  disturbance.  Calculations 
performed  for  models  similar  to  that  in  Fig.  lc  but  with  more  diffuse 
boundaries  (larger  Ar)  produce  results  similar  to  those  in  Fig.  8 except 
for  the  diminished  importance  of  reflections  and  the  corresponding  reduc- 
tion of  the  standing  wave  Intensity. 

A final  purpose  of  this  analysis  is  to  assess  whether  lateral  nonuni- 
formities Impair  the  accuracy  of  WKB  calculations  of  the  effects  of  high-alti- 
tude nuclear  detonations  on  ELF  propagation.  Field  (.1978)  evaluates  the 
effects  of  bursts  at  midpath;  here  we  focus  on  the  effect  of  a burst 
detonated  over  the  transmitter.  Bursts  at  altitudes  below,  say,  150  km 
produce  ionospheric  disturbances  with  lateral  scales  up  to  perhaps  1500  km 
(e.g..  Field  and  Engel,  1965).  The  analytic  models  used  to  calculate 
Figs.  2 through  4 are  reasonably  representative  of  that  range  of  parameters. 

To  supplement  those  results,  and  to  examine  less  idealized  nuclear 
environments.  Fig.  9 shows  calculated  results  for  two  nominal  bursts 
producing  1)  two  megatons  of  fisison  debris  at  300  km  altitude  1 minute 
after  detonation,  and  2)  two  megatons  of  fission  debris  at  1000  km 
altitude  10  minutes  after  detonation.  The  values  of  S(r)  used  as  inputs 
to  Eq.  (13)  were  calculated  by  Greifinger  and  Greifinger  (1977).  Even  for 
the  low  assumed  frequency  of  45  Hz,  Fig.  9 indicates  that  the  WKB  solution 
affords  an  excellent  approximation  of  the  full-wave  solution. 


-Amplitude  of  W versus  distance  for  propagation  from  undisturbed  to 
disturbed  uniform  regions  separated  by  diffuse  boundary  at  r = 5 N 


Amplitude  of  W versus  distance  for  two  nominal  nuclear  disturbances  over 
transmitter:  daytime,  f * 45  Hz. 
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IV.  CONCLUSIONS 

Recasting  the  equation  for  the  lateral  dependence  of  fields  into  an 
integral  equation  is  a convenient  means  of  obtaining  numerical  solutions 
for  long-wave  propagation  under  nonstratif led  conditions.  The  integral 
version  accounts  for  most  full-wave  properties,  including  diffraction 
around  a localized  disturbance  and  reflection  from  lateral  gradients. 
Numerical  solutions  based  on  model  disturbances  having  lateral  gradients 
in  the  direction  of  propagation  reveal  a standing  wave  pattern  in  front 
of  the  disturbance.  The  pattern  is  pronounced  if  the  waveguide  properties 
change  substantially  over  a distance  equal  to  about  one-sixth  of  a wave- 
length, and  is  minor  if  the  disturbance  is  more  diffuse. 

The  often  used  WKB  approximation  of  course  ignores  diffraction  and 
gradient  reflection.  Thus  it  omits  the  standing  wave  pattern,  thereby 
giving  poor  accuracy  for  regions  in  front  of  a disturbance  having  a 
relatively  abrupt  boundary.  Because  it  ignores  gradient  reflections, 
the  WKB  approximation  also  overestimates  the  signal  transmitted  beyond 
a nonuniformity  in  the  waveguide.  However,  for  all  models  considered — 
abrupt  or  diffuse — the  overestimate  is  minor. 

The  integral-equation  method  developed  here  is  valid  at  any  frequency 
for  which  waveguide  modes  describe  terrestrial  propagation.  Nonetheless, 
its  practical  utility  is  probably  limited  to  ELF  since  the  lateral  proper- 
ties of  the  earth-ionosphere  waveguide  are  usually  gradual  enough  to 
permit  use  of  the  WKB  approximation  at  higher  frequencies. 
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Appendix  A 

INTEGRAL  EQUATION  FOR  AN  AZ IMUTHALLY  SYMMETRIC  DISTURBANCE 


If  the  disturbance  is  cylindrically  symmetric  and  centered  at  the 

2 
0 


2 2 

origin,  S (r)  - S.  ■ F(r),  and  Eq.  (9)  (p.  7)  can  be  rewritten 


2 °°  2ir 

iKr,0)  - «|»o(r,0)  - J"  r'dr'F(r')  \ dO'Hg 


(2) 


(kSor2H(r\0')  , (A-l) 


where  r and  0 are  the  usual  circular  cylindrical  coordinates,  and 


r2  - (r')2  + r2  - 2rr'  cos(0'  - 0) 


(A-2) 


To  obtain  a one-dimensional  integral  equation  in  the  variable  r,  we  expand 
the  terms  in  Eq.  (A-l)  in  powers  of  cos0  as  follows: 


«Kr 


,9> " S c"*" 


cos  m0 


(A-3) 


OO 

H02>(kS0r2)  "2^  en}Jn(kSOr<)  Hn2)(kS0r>)JCO8n(9' "0)  • (A"A) 


,(2) 


*ft(r,0)  - SnHv  '(kSftr)  cos  6 


(A-5) 


where  ■ 1 or  2 for  n ■ 0 or  n > 1,  respectively,  and  r<  and  r>  denote, 

C respectively,  the  lesser  or  greater  of  r and  r'.  Equation  (A-4)  is  the 

addition  theorem  for  Hankel  functions  (e.g. , Magnus  and  Oberhsttinger , 1943), 
and  Eq.  (A-S)  is  the  well-known  solution  for  a laterally  homogeneous, 

G undisturbed  ionosphere  (see  Eq.  (3),  p.  4). 
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By  inserting  Eqs.  (A-3)  through  (A-5)  into  Eq.  (A-l),  and  noting  that 


2tt  /n  cos0  If  n ■ n 

(dB'cosnO'cosmfQ'  - 8)  ■ \ , 

( 0 if  n ft  m 


(A-6) 


it  follows  that 


QO 


cosn»0  » SqH^2^  (kS^r)  cos0  - cos0 


00  00 
nrO  m J 


r'dr'F(r') 


,(2) 


* VkV<>  m (kSOr>)  ^ * 


(A-7) 


Using  trigonometric  orthogonality  relations,  it  can  be  shown  that  only 
is  nonzero,  and  the  equation  for  the  radial  dependence  of  <|>  becomes 


>j(r)  - Sqh|2^ (kSQr) 


ink 

2 


00 

J r'dr'F(r' 


) ^kSor<)Hj2)  (kSQr>)i('1  (r')  . (A-8) 


The  relative  propagation  function,  W,  is  defined  by 

Mr)  - <|»Q(r)W(r) 


(A-9) 


and  combining  Eqs.  (A-3),  (A-5),  (A-8),  and  (A-9)  gives  the  following 
result: 


W(r)  - 1 - 


r'dr'F(r’) 


(kS0r<)H{2)  (kSQr>)H{2)  (kS^* ) 


H 


(2) 


(kSQr) 


W(r') 


(A-10) 


Aside  from  some  differences  in  notation,  Eq.  (A-10)  is  identical  to 
Eq.  (13)  (p.  9). 
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Appendix  B 

WKB  SOLUTION  FOR  AN  AZIMUTHALLY  SYMMETRIC  DISTURBANCE 


If  S depends  only  on  the  distance,  r,  from  the  origin,  Eq.  (5)  can  be 


written 


"S'  ♦ + r ^ + (k2S2(r)  - 1/r2)^  - 0 

dr  1 


By  making  substitution  ip  ■ u/r  , Eq.  (B-l)  can  be  transformed  into 


the  following  standard  form: 


sV)  --4 

dr  L 4k z 


If  S >>  1/kr,  Eq.  (B-2)  is  identical  to  the  equation  of  a plane  wave 
propagating  in  a medium  of  refractive  index  S,  and  the  WKB  solution 
normalized  to  the  definition  of  ip  given  in  Eqs.  (3)  and  (4)  can  be 


written  (e.g. , Budden,  1961b): 


, *s  S(r  ■ 0) 

<1*  M u/r — q—  exp 

r S(r) 


h / 


S(r')dr' 


Using  the  definition  W * , the  WKB  approximation  to  the  relative 


propagation  function  is  found  to  be 


exph  /is<r,>  - s°H  • 


which  is  Eq.  (14)  (p.  9). 


(B-l) 


(B-2) 


( B— 3) 


(B-4) 
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